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Abstract 

The processes e~'~e~ ti/hh+Wiggs allow to measure the Yukawa couplings between 
Higgs bosons and heavy quarks in supersymmetric theories. The complete set of 
next-to-leading order SUSY-QCD corrections to the cross sections of these processes 
have been determined in the minimal supersymmetric extension of the Standard 
Model. They turn out to be 0(10 — 20%) and thus important for future linear e^e~ 
colliders. 
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1 Introduction 



The Higgs mechanism is a cornerstone of the Standard Model (SM) and its supersym- 
metric extensions jT]. In the minimal supersymmetric extension of the SM (MSSM) two 
isospin Higgs doublets have to be introduced in order to generate masses of up- and 
down-type fermions After electroweak symmetry breaking three of the eight degrees 
of freedom are absorbed by the Z and W gauge bosons, implying the existence of five 
elementary Higgs particles. These consist of two CP-even neutral (scalar) particles 
one CP-odd neutral (pseudoscalar) particle A and two charged bosons H^. At leading 
order the MSSM Higgs sector is fixed by two independent input parameters which are 
usually chosen to be the pseudoscalar Higgs mass Ma and tg/3 = v^jvx, the ratio of the 
two vacuum expectation values. Including the one-loop and dominant two-loop correc- 
tions the upper bound of the light scalar Higgs mass is < 135 GeV 3j. The couplings 
of the various neutral Higgs bosons to fermions and gauge bosons, normalized to the SM 
Higgs couplings, are listed in Table ^ where the angle a denotes the mixing angle of the 
scalar Higgs bosons h, H. An important property of the bottom Yukawa couplings is their 
enhancement for large values of tg(3, while the top Yukawa couplings are suppressed for 
large tgjS. This implies that Higgs radiation off bottom quarks plays a significant role 
at linear e~^e~ colliders in the large tg/3 regime, while Higgs radiation off top quarks is 
relevant for small and moderate values of tg/?. 
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Table 1: Higgs couplings in the MSSM to fermions and gauge bosons [V = W, Z] relative 
to the SM couplings. 

For the computation of the SUSY-QCD corrections we need the Higgs couplings to 
stop and sbottom squarks in addition, as well as the squark couplings to the photon and 
Z boson. The scalar superpartners fi^R of the left- and right-handed fermion components 
mix with each other. The mass eigenstates /i^2 of the sfermions are related to the current 
eigenstates /^^r by mixing angles 9f, 

fi = fLCOsOf + fRsinOf 

h = -fLsinOf + fRCOsOf . (1) 
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The mass matrix of the sfermions in the left-right-basis is given by |3]^ 

+ mj mf{Af - firj) 



(2) 



with the parameters = l/r„ = tg/3. The parameters Aj denote the trihnear scalar 
couplings of the soft supersymmetry breaking part of the Lagrangian. The mixing angles 
acquire the form 

2mf(Af-arf) M} - M? 

/I 72 Jl J2 

and the masses of the squark mass eigenstates are given by 

1 r 



Since the mixing angles are proportional to the masses of the ordinary fermions, mixing 
effects are only important for the third-generation sfermions. The neutral Higgs couplings 
to sfermions read |5| 

g'f^f^ = mjgf + MlefSin'^ewgt 

4;, = -^{f^gt - Afgt) , (5) 

with the couplings gf listed in Table El hf denotes the third component of the electroweak 
isospin, Cf the electric charge of the fermion / and sin 6w the Weinberg angle, rrif the 
fermion mass and Mz the Z-boson mass. 

Finally the sfermion couplings to photons are given by j3] 



and to Z bosons by 



= ^^f^^^ (6) 



gf - = ihf-ef sin^ Ow) 



IrSr 



SwCw 



w 



4/. = 0- (7) 

All these couplings have to be rotated to the mass eigenstates by the mixing angle 6f. 
The coupling e = y/Ana denotes the elementary electric charge. The Yukawa couplings 

^For simplicity, the _D-terms have been absorbed in the sfermion-mass parameters M? 
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Table 2: Coefficients of the neutral MSSM Higgs couplings to sfermion pairs. 



of squarks to gluinos and quarks are uniquely determined by the strong coupling as- The 
corresponding Feynman rules can be found in Ref. jlj. 

Neutral Higgs radiation off top or bottom quarks [Q = t, b] in e~^e~ collisions, 

e+e-^gg/ = h,H, A], (8) 

is a suitable process for measuring the Yukawa couplings in supersymmetric theories jB], 
particularly for the light Higgs boson h and for moderately heavy Higgs bosons H and 
A. In the following we present the cross sections for these processes including the next- 
to-leading order (NLO) SUSY-QCD corrections. The pure QCD corrections have been 
determined before in the SM j7] and the MSSM ^ (approximate results can be found in 
Refs. |9llUj ). They significantly modify the cross section at high collider energies. For SM 
Higgs boson radiation off top quarks the full electroweak corrections have recently become 
available jTTj. They turn out to be of similar magnitude as the pure QCD corrections so 
that an analogous result may be expected for MSSM Higgs bosons. First results of the 
SUSY-QCD corrections to tth production appeared in Ref |12j. We will comment on the 
large size of these corrections in this paper. 

2 Calculation of NLO SUSY-QCD Corrections 

In this section we will set up our notation and summarize the present status of the pro- 
cesses under investigation, before we will describe the novel SUSY-QCD corrections in 
detail. 

2.1 Leading-Order Cross Sections and QCD Corrections 

At leading order (LO) Higgs radiation off heavy quarks is described by the Feynman di- 
agrams of Fig. n It splits into three different classes of contributions: (i) Higgs radiation 
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Figure 1: Individual channels of scalar and pseudoscalar Higgs radiation off heavy quarks 
Q = t,b in e^e^ collisions. 

off the heavy (anti) quark, (ii) Higgs radiation off the Z boson (only scalar Higgs radi- 
ation) and (iii) Z boson splitting into scalar-pseudoscalar Higgs pairs with one of them 
dissociating into a heavy QQ pair. Depending on the masses of the corresponding par- 
ticles resonant contributions will arise, which require the inclusion of finite decay widths 
of the Z and Higgs bosons in the corresponding propagators. We have used conventional 
Breit-Wigner propagators for the resonant Z ^ bb and 0° — * ti/bb decays as in previous 
analyses 8j. The LO matrix elements can be decomposed according to the different Higgs 
couplings, 

Mto = g^'D^'' + g'g'D^'' + gS-^^"" . (9) 

The LO cross sections can be cast into the form 

aLo = J dPSs E \Mto\' (10) 

spins,colours 

with dPSs denoting the three-particle phase space element. The sum has to be taken over 
the spins of the initial and final-state fermions and colours supplemented by an average 
over the initial e~^e~ spins. 

We include the QCD corrections of Ref. [S] with the QCD coupling evaluated at 
NLO with 5 active flavours at the renormalization scale hr = ^/s with s being the e~^e~ 
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cm. energy squared. The coupling is normalized to a;s(M|) = 0.119. The bottom 
Yukawa couplings are computed at the scale of the corresponding Higgs-momentum flow. 
This choice absorbs large logarithmic contributions of the pure QCD corrections |8jj. 

The value of the electromagnetic coupling is taken to be a = 1/128 and the Weinberg 
angle to be = 0.23. The mass of the Z boson is set to Mz = 91.187 GeV, and the pole 
masses of the top and bottom quarks are set to rrit = 175 GeV^ and rrib = 4.62 GeV^, 
respectively [S]. The masses of the MSSM Higgs bosons and their couplings are related 
to tgP and the pseudoscalar Higgs-boson mass Ma- In the relations we use, higher-order 
corrections up to two loops in the effective-potential approach are included |T3I. The Z 
boson width has been chosen as Tz = 2.49 GeV, and the Higgs boson widths have been 
computed with the program HDECAY |15j . 

2.2 SUSY-QCD Corrections 

The NLO SUSY-QCD corrections arise from virtual gluino and stop/sbottom exchange 
contributions as depicted in Fig. El They consist of self-energy, vertex and box contribu- 




Figure 2: Typical diagrams of the NLO SUSY-QCD corrections to e~^e QQ(fP [Q = 
t, h] mediated by gluino g and squark Q = t,b exchange. 



^The top mass has been chosen in accordance with the definitions of the Snowmass-benchmark points 
of the MSSM |13|. 

■^This value for the perturbative pole mass of the bottom quark corresponds in NLO to an MS mass 
mbirrib) = 4.28 GeV. 
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tions, which are calculated within dimensional regularization in the standard way. Since 
all virtual particles are massive, no infrared nor collinear singularities arise. The ultra- 
violet divergences are removed by the renormalization of the quark masses and Yukawa 
couplings. The latter is connected to the quark mass renormalization. In the case of 
tt+Higgs production the top mass has been renormalized on-shell in the propagators as 
well as in the Yukawa couplings. The same prescription has also been chosen for the 
bottom mass, since the virtual gluino and sbottom masses are too large to develop large 
logarithmic contributions. Thus the renormalization of the bottom quark mass is given 
by 

f^b = [1 + Sqcd + Ssqcd] 

Sqcd = -^Tie)iAnY 

TT 

S(mb) 

OSQCD - _ / ' 

where denotes the bare bottom mass, Tni,{jj^) the MS mass at the scale and S(^s)qc'd 
the corresponding (SUSY-)QCD counter terms. The SUSY-QCD contribution to the 
bottom quark self-energy is displayed by T,{mb) with on-shell momentum. This renor- 
malization prescription ensures that the gluino and sbottom contributions are decoupled 
from the running of the bottom Yukawa couplings. Thus we are left with the pure MS 
Yukawa couplings of QCD. 

The final result can be cast into the form 

a(e+e- ^ QQcj)') = aLo{e+e- ^ QQ/) |l + [Cqcd + Csqcd] ^| , (12) 

where aLo{.^^^~ QQ4'^) denotes the LO cross section and C(^s)qcd the coefficients of 
the (SUSY-) QCD corrections. 

For large values of tg/3 there are significant non-decoupling corrections to hhcf)^ pro- 
duction, which can be absorbed in the bottom Yukawa couplings in a universal way |16j^. 
In Refs. jl8ll9j it has been shown that these contributions can be resummed to improve 
the reliability of the perturbative result. The resummed bottom Yukawa couplings are 
given by^ 

~9'b = 
= 

9t = T^(l-£i^) (13) 




"'it should be noted that these contributions vanish for large sbottom and gluino masses while keeping 
the fi parameter fixed 03 . Non-decoupling effects only arise, if the fi parameter is increased together 
with the SUSY particle masses. 

^Analogous effective couplings can be defined for top quarks, too, but in this case the non-decoupling 
contributions are small and thus do not require resummation. 



6 



with the quantities [121 
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(a — b){b — c)(c — a) 

If the LO cross sections are expressed in terms of these resummed bottom Yukawa cou- 
phngs, the corresponding NLO pieces have to be subtracted from the coefficient Csqcd to 
avoid double counting. This is equivalent to an additional (finite) renormalization, given 
explicitly by 
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(15) 



Thus we have to add additional finite counter terms to the virtual matrix elements^ 

_ h/ 
'■SQCD — 9b 



AM'l'^rr. = .'J^'ArC^+gtAtC 



A-^5QCD = 9t^t'^i+9l^lV^ + 9b^b'^z. 
In the LO matrix elements we use the resummed Yukawa couplings, 



-hlHpLO I h/HpLO , nApLO 
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(16) 



(17) 



so that the SUSY-QCD corrections to the cross sections are given by 

a. 



AcTsQCD — (^LoCsQCD- 



dPS, 2^e Y: Mt^MfQCD 

spins , colour s 



In the QCD corrections we also insert the resummed bottom Yukawa couplings every- 
where, since the non-decoupling terms A^ factorize from the pure QCD corrections ini- 
tiated by light particle interactions. 

^Note that in the matrix elements of the SUSY-QCD corrections we keep the unresummed bottom 
Yukawa couphngs in order to avoid artificial singularities for vanishing mixing angle a |19|. 



3 Results 



The numerical results will be presented for a linear e~^e~ collider with cm. energy of 
1 TeV. We have chosen the Snowmass point SPSS for Higgs radiation off top quarks and 
SPSlb for the bottom quark case The MSSM parameters of these two benchmark 
scenarios are given by^ 

SPSS: 



tgf3 = 
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The pseudoscalar Higgs mass is left free in both scenarios in order to scan the correspond- 
ing Higgs mass ranges. 

The total cross section for pseudoscalar Higgs radiation off top quarks is displayed at 
LO and NLO in Fig. ISt- The cross section is small for pseudoscalar Higgs masses below 
about 3S0 GeV, while above it rapidly increases to a level of 1 fb due to the intermediate 
on-shell H ti decay. The total size of the corrections amounts to about (9(10%) apart 
from the threshold of the resonant contribution, where the Coulomb singularity raises 
the QCD corrections to more than 100% |i7j. The Coulomb singularity is an artefact of 
the narrow-width approximation. A proper treatment of the threshold region requires 
the inclusion of finite-width effects and QCD-potential contributions. Thus the result 

^We have neglected the corresponding translations of DR masses into MS masses, since they are not 
relevant for the characterization of the results. 
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obtained in this work is not valid in a small margin around the tt threshold of the resonant 
part. The individual relative corrections, defined as aNio = c"lo(1 + ^qcd + ^sqcd), can 
be inferred from Fig. ^jp. Except for the threshold region of the resonant part, the QCD 
corrections are of moderate size 0. The SUSY-QCD corrections are of similar magnitude 
as the pure QCD corrections but of opposite sign. Thus, we observe a large cancellation 
of the QCD corrections against the SUSY-QCD part. This signalizes the importance of 
including both types of corrections in future analyses. An analogous picture emerges for 
the light and heavy scalar Higgs bosons as shown in Figs. and b. The cross section 
for the light scalar Higgs boson is always of the order of 1 fb with small corrections 
due to the partial cancellation of QCD and SUSY-QCD corrections (see Fig. HJd). The 
QCD Coulomb singularity for ~ 350 GeV is much more pronounced than in the 
pseudoscalar case, since for the heavy scalar Higgs boson the S-wave pseudoscalar Higgs 
decay A ^ tt constitutes the resonant part^. The relative threshold corrections remain 
finite in both cases due to the remaining continuum contributions. 

The results for bbA production are presented in Fig. El The total cross section, shown in 
Fig.EK; reaches a size of (9(10 fb) for smaller pseudoscalar masses. The relative corrections 
are depicted in Fig. The pure SUSY-QCD and total corrections are shown without 
and with resummation of the A;, terms according to Eqs. (jl3flTH|) . It is clearly visible that 
the resummed bottom Yukawa couplings absorb the bulk of the SUSY-QCD corrections, 
so that the terms of Eq. (fT^ provide a reasonable approximation of the final result. After 
resummation the SUSY-QCD corrections cancel against the pure QCD corrections to a 
large extent. Thus, as in the top quark case the inclusion of both corrections is of vital 
importance. A comparison of the total resummed and unresummed NLO cross sections 
in Fig. implies good agreement within 10% and thus a significant improvement of the 
perturbative stability from LO to NLO. An analogous picture emerges for the light and 
heavy scalar Higgs bosons as can be inferred from Fig. IHl Again the resummed Yukawa 
couplings absorb the bulk of the SUSY-QCD corrections. A significant cancellation of the 
QCD and SUSY-QCD corrections is observed after resummation in this case, too. The 
drops of the relative corrections towards ~ 500 GeV in Figs. El and El are caused by 
the kinematical closure of the intermediate on-shell HA pair production. 

Higgs radiation off bottom quarks is, however, dominated by the resonant h,H bb 
decays in the pseudoscalar case and the resonant Z,A-^bb decays in the scalar case. 
Thus, the absorption of the bulk of the SUSY-QCD part by the resummed Yukawa cou- 
plings could be expected from the analogous findings for the corresponding Higgs decays 
[T^ . In order to investigate, if this also holds for continuum bb(j)^ production, we analyze 
the Higgs energy distribution in Fig. [71 for pseudoscalar Higgs radiation off bottom quarks 
for a pseudoscalar Higgs mass Ma = 200 GeV. The dimensionless parameters x^o are 
defined as x^o = 2E^o/y^. Fig. |2k displays the xa distribution at LO and NLO, while 
Fig. 03 exhibits the individual relative corrections. The sharp peak at a:^ ~ 1 originates 
from the resonant h, H ^ bb decays, while the regions apart from the peak represent 

^ Since the pseudoscalar Higgs boson A carries the same quantum numbers as the ground state 
of the tt pair, the decay A tt is dominated by an S'-wave contribution at threshold. In contrast the 
scalar Higgs decay H ^ tt suffers from a P-wave suppression at threshold. 
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Figure 3: (a) Production cross sections of pseudoscalar Higgs radiation off top quarks in 
e'^e~ collisions. The LO cross section is depicted by the dashed line and the full QCD- and 

SUSY-QCD corrected cross section by the full line; (b) Relative QCD, SUSY-QCD and 
total corrections to pseudoscalar Higgs radiation off top quarks. The sharp (finite) peak 
around Ma = 350 GeV originates from the Coulomb singularity in the QCD corrections 
to the resonant H ^ tt decay. 
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Figure 4: (a) Production cross sections of heavy and light scalar Higgs radiation off top 
quarks in e'^e" collisions. The LO cross section is depicted by the dashed line and the full 
QCD- and SUSY-QCD corrected cross section by the full line; (b) Relative QCD, SUSY- 
QCD and total corrections to scalar Higgs radiation off' top quarks. The sharp (finite) peak 
around Mh = 350 GeV originates from the Coulomb singularity in the QCD corrections 
to the resonant A ^ tt decay. 
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Figure 5: (a) Production cross sections of pseudoscalar Higgs radiation off bottom quarks 
in e+e~ collisions with (red curves) and without (black curves) resummation of the Af, 
terms. The LO cross section is depicted by the dashed line and the full QCD- and SUSY- 
QCD corrected cross section by the full line; (b) Relative QCD, SUSY-QCD and total 
corrections to pseudoscalar Higgs radiation off bottom quarks with (red lines) and without 
(black lines) resummation. The pure QCD corrections are indistinguishable in both cases. 
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Figure 6: (a) Production cross sections of heavy and light scalar Higgs radiation off bottom 
quarks in e~^e~ collisions with (red curves) and without (black curves) resummation of the 
Afe terms. The LO cross section is depicted by the dashed line and the full QCD- and 
SUSY-QCD corrected cross section by the full line; (b) Relative QCD, SUSY-QCD and 
total corrections to scalar Higgs radiation off bottom quarks with (red lines) and without 
(black lines) resummation. The pure QCD corrections are indistinguishable in both cases. 
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continuum bbA production. The resulting picture indeed turns out to be analogous to 
the total cross sections. The bulk of the SUSY-QCD corrections can be absorbed by the 
resummed bottom Yukawa couplings leaving moderate residual corrections. These cancel 
the pure QCD corrections to a large extent in the resonant as well as the continuum 
regions. 

The light scalar Higgs energy distribution for Higgs radiation off top quarks is shown 
in Fig. IHl for a light scalar Higgs mass = 100 GeV. For Xh ^ 0.8 both the QCD 
and SUSY-QCD corrections are of moderate size and cancel each other partly. The sum 
of the corrections amounts to a few per cent. The sharp rise of the QCD corrections 
towards Xh ~ 0.9 is induced by the Coulomb singularity at the subthreshold of the tt 
pair S . It leads to a finite cross section at the upper bound of the Xh range. Since the 
total corrections are not constant, the shape of the Higgs energy distribution is slightly 
modified from LO to NLO as can be inferred from Fig. |Hk- 

We performed a comparison of our tih results with those of Ref. [12]. We could not 
tune our MSSM scenario exactly to the one of Ref. since we do not use the same 
MSSM routines to calculate the SUSY-particle masses and couplings. Nevertheless, we 
have found reasonable agreement with the results in a scenario very close to the one in 
Ref. ^2] ■ However, this scenario is characterized by a very extreme choice of parameters. 
The light stop ti is very light (m^^ ~ 100 GeV), and the gluino mass of 200 GeV is chosen 
such that rrit + Mh ~ rrig + m^^, so that the thresholds of the virtual Yukawa vertex 
corrections in the momentum of the off-shell (anti)top quarks, which split into t/t + h, are 
very close to the average momentum flow through the (anti)top propagator (see the third 
diagram of Fig. This enhances the size of the SUSY-QCD corrections considerably. 
However, the theoretical uncertainties in the threshold regions are large, and the results 
cannot be trusted quantitatively. A better perturbative treatment requires the inclusion 
of finite-width effects of the unstable particles as well as QCD-potential effects. Since the 
scenario of Ref. J2] belongs to a very particular region in the MSSM parameter space, 
this result does not represent the overall size of the corrections. It can clearly be inferred 
from Fig. 6 of Ref. jl2j that the corrections are of moderate size away from this threshold 
region. 

4 Conclusions 

We have presented the full SUSY-QCD corrections to neutral MSSM Higgs radiation 
off top and bottom quarks at linear e'^e~ colliders. The size of the corrections is of 
0(10 — 20%) and of similar magnitude as the pure QCD corrections obtained in the past, 
but of opposite sign, so that significant cancellations occur. This underlines the relevance 
of including these corrections in future analyses of these processes at linear e~^e~ colliders. 

At large values of tg/3 Higgs radiation off bottom quarks provides a possibility to 
measure tg(3 [201 ■ In the past it has been demonstrated that the bulk of the pure QCD 
corrections can be absorbed in the running bottom Yukawa couplings, defined at the 
scale of the corresponding Higgs momentum flows jS]. The SUSY-QCD corrections on 
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Figure 7: (a) Higgs energy distribution of pseudoscalar Higgs radiation off bottom quarks 
in e~^e~ collisions with (red curves) and without (black curves) resummation of the At 
terms. The LO cross section is depicted by the dashed line and the full QCD- and SUSY- 
QCD corrected cross section by the full line. The peak at ^ 1 originates from the 
resonant h,H ^ bb decays; (b) Relative QCD, SUSY-QCD and total corrections to pseu- 
doscalar Higgs radiation off bottom quarks with (red lines) and without (black lines) re- 
summation. The pure QCD corrections are indistinguishable in both cases. 
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Figure 8: (a) Higgs energy distribution of light scalar Higgs radiation off top quarks in 
e'^e~ collisions. The LO cross section is depicted by the dashed line and the full QCD- 
and SUSY-QCD corrected cross section by the full line; (b) Relative QCD, SUSY-QCD 
and total corrections to light scalar Higgs radiation off top quarks. The sharp rise of the 
QCD corrections towards Xh ~ 0.9 is induced by the Coulomb singularity. 
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the other hand are dominated by the non-decouphng Am;, terms, which can be absorbed 
and resummed in the corresponding bottom Yukawa couphngs. We have shown that this 
absorption reduces the SUSY-QCD corrections to a moderate size, but still cancels the 
pure QCD corrections to a large extent in the resonant as well as continuum regimes. 

Based on the size of the electroweak corrections to Standard Model UH production 
the determination of the full SUSY-electroweak corrections to these processes turns out 
to be necessary, to reduce the overall theoretical uncertainties to a level which allows 
accurate measurements of the Higgs Yukawa couplings. 
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